Raman scattering study of delafossite magnetoelectric multiferroic 

compounds: 
CuFe02 and CuCrOa 

Oktay Aktas, Kim Doan Truong, Tsuyoshi Otani, Geetha Balakrishnan, 
Maynard J. Clouter, Tsuyoshi Kimura, and Guy Quirion 

October 5, 2011 

Abstract 

Ultrasonic velocity measurements on the magnetoelectric multiferroic compound CuFe02 reveal that 
the antiferromagnetic transition observed at Tjvi = 14 K might be induced by an R3m C2/m pseudo- 
proper ferroelastic transition |1J. In that case, the group theory states that the order parameter associated 
with the structural transition must belong to a two dimensional irreducible representation Eg {x^ — , 
xy). Since this type of transition can be driven by a Raman Eg mode, we performed Raman scattering 
measurements on CuFe02 between 5 K and 290 K. Considering that the isostructural multiferroic com- 
pound CuCr02 might show similar structural deformations at the antiferromagnetic transition Tjvi = 
24.3 K, Raman measurements have also been performed for comparison. At ambient temperature, the 
Raman modes in CuFe02 are observed at oje, = 352 cm"^ and cjAg ~ 692 cm^^, while these modes 
are detected at ujEg ~ 457 cm~^ and uj Ag = 709 cm~^ in CuCr02. The analysis of the temperature 
dependence of modes shows that the frequency of all modes increases down to 5 K. This typical behavior 
can be attributed to anharmonic phonon-phonon interactions. These results clearly indicate that none 
of the Raman active modes observed in CuFe02 and CuCr02 drive the pseudoproper ferroelastic tran- 
sition observed at the Neel temperature Tjvi. Finally, a broad band at about 550 cm~^ observed in the 
magnetoelectric phase of CuCr02 below Tn2 could be attributed to a magnon mode. 

1 Introduction 

CuFe02 and CuCr02 belong to the delafossite frustrated antiferromagnets with the chemical formula ABO2 
in which A is a nonmagnetic monovalent ion (Cu and Ag) while B is a magnetic trivalent ion such as Fe and 
Cr [2 [3l m [5] . Some of these compounds, including AgCr02, CuFe02, and CuCr02, belong to the trigonal 
R3m space group at room temperature and undergo a series of magnetic phase transitions [3l lU |5] at low 
temperatures as a result of geometrical frustration of magnetic ions sitting on a triangular lattice. 

In the case of CuFe02, two antiferromagnetic transitions are observed at zero field. In its ground state, 
Fe"'"'' ions order into a collinear commensurate four-sublattice (ttii) structure, while between = 11 K 
and T/vi = 14 K, the magnetic order is incommensurate with the magnetic moments also pointing along the 
c-axis [6]. With the application of a field parallel to the c-axis, a series of new magnetic orders is stabilized 
below Tn2- Between 7 T and 13 T, CuFe02 shows a proper screw spin configuration where the spins lie in 
the R3m mirror plane perpendicular to the magnetic modulation vector q|| [110] (hexagonal basis) (HHHTIIH]. 
At higher fields, several other spin configurations are observed: a c-axis collinear 5-sublattice (ttt4-4-) state 
(13 T < H < 20 T), a c-axis collinear 3-sublattice (fU) structure (20 T < H < 34 T), a canted 3-sublattice 
state (34 T < H < 49 T), and a noncolinear incommensurate spin- flop phase which is close to the 120° spin 
structure for 49 T < H < 70 T, followed by a transition to the paramagnetic state at 70 T. [U [9] . 

While CuCr02 is isostructural to CuFe02 at room temperature, its magnetic phase diagram is signifi- 
cantly different [31 [51 [TDj . According to specific heat and magnetic susceptibility measurements ^ , CuCr02 
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shows anomalies at Tjvi = 24.3 K and T7V2 — 23.6 K. The magnetic order in the intermediate temperature 
range T/vi < T < T]si2 is interpreted as a coUinear state with S||c 0, while recent neutron diffraction mea- 
surements 111 ( 112 ) reveal an incommensurate proper screw spin structure with q||[110] below Ti^2- This spin 
configuration is very similar to the one observed in CuFe02 between 7 T and 13 T. Moreover, additional 
studies on both compounds |U [31 [S] show that an electric polarization P|j[110] is only induced upon the 
emergence of this proper-screw spin order. Under this scenario, the usual inverse Dzyaloshinskii-Moriya 
(DM) interaction P ~ Vij x (Si x Sj) [131 Ej cannot account for the induced polarization as the q-vector 
of the spin modulation is perpendicular to the spiral-plane. An alternative possibility, proposed by Arima 
et al. |15j . is that the polarization is induced by spin-orbit coupling. Thus, CuFe02 and CuCr02 represent 
a different class of magnetoelectric multifcrroics in which the mechanism leading to the magnetoelectric 
coupling is still uncertain. 

Other particular properties of CuFe02 have also been recently revealed via sound velocity measure- 
ments [1] m [16] . These measurements show softening on specific elastic constants as the temperature is 
decreased down to T/vi = 14 K. The data analysis indicates that this peculiar behavior is characteristic of 
an _R3m ^ C2/m pseudoproper ferroelastic transition, consistent with neutron |17j and x-ray diffraction 
measurements. Furthermore, according to the group theory [TH], the order parameter associated with the 
structural transition must belong to a two dimensional irreducible representation (IR) Eg {x^ — , xy). 
As none of the spin components belong to this IR, these measurements indicate that the magnetic order 
in CuFe02 is stabilized by the ferroelastic structural transition, putting in evidence the role played by the 
spin-lattice coupling in this family of multiferroic materials. Thus, the true origin of the structural transition 
observed at T/vi remains a mystery. One possibility is that the transition is driven by a Raman mode as in 
other pseudoproper ferroelastic materials [THl [201 [HI [12] • Regarding isostructural CuCr02 , recent megne- 
tostriction measurements |23j show evidence for a structural phase transition at Tjvi = 24.3 K. Furthermore, 
as preliminary sound velocity measurements on CuCr02 (shown in Fig. [ij reveal softening similar to that 
observed in CuFe02 [1|, the transition observed at T^i — 24.3 K in CuCr02 might also be ferroelastic. In 
order to possibly identify the order parameter associated with these pseudoproper ferroelastic transitions, 
we performed Raman scattering measurements on the isostructural compounds CuFe02 and CuCr02. The 
remainder of the paper is organized as follows. Experimental methods are discussed in Sec. [2] while results 
and discussion are presented in Sec. [3] Finally, conclusions are made in Sec. |4] 

2 Experiment 

Single crystals of CuFe02 were grown by the floating zone method using a four mirror image furnace [24) . 
CuFe02 samples used in the measurements had an area of ^^2 mm x ~2 mm and was ~1 mm long along 
the c-axis. Single crystals of CuCr02 were grown from Bi203 flux [5]. The samples were platelets with a 
length of 0.4 mm along the c-axis. The surface area was approximately 2 mm x 2 mm. Prior to Raman 
scattering experiments, samples were polished using abrasive slurry with 50 nm AI2O3 grains in order to 
minimize surface scattering. Room temperature Raman measurements were performed using two different 
experimental setups. Using an Ar+ laser operating at 514.5 nm, the Raman spectra were collected by a 
double grating spectrometer (Spex Industries, model 1401), a photomultiplier tube (Perkin Elmer, MP 900 
series), and a photon counter (Princeton Applied Research, model 1109). For CuCr02, 28 mW of exciting 
beam power was used while it was increased to 50 mW to obtain the spectra of CuFe02. For cross polarization 
measurements on CuFe02, the beam power was increased up to 100 mW which caused local heating on the 
sample. In the second setup, the 0.5 cm"'^ resolution micro-Raman measurements were performed with a 
632 nm He-Ne laser, a double grating spectrometer (Jobin Yvon, model Labram-800) and a liquid-nitrogen 
cooled CCD detector. In order to minimize sample heating, 0.3 mW incident beam power with a 3 /zm spot 
size (4000 W/cm^) was used. Mode parameters were obtained by a fit to the data using Lorentzian functions 
for the observed modes. 
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Figure 1: Relative variation of the elastic constant Cqq in CuCr02 as a function of temperature obtained by 
sound velocity measurements. Cee shows a 30 % reduction at Tjvi = 24.3 K relative to the value at 150 K, 
indicating that the antiferromagnetic transition at Tjvi might also be ferroelastic. 
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3 Experimental results and discussion 



3.1 Raman spectra of CuFe02 and CuCr02 at room temperature 

Delafossite compounds (space group i?3m) such as CuFe02 and CuCr02 have one formula unit per unit 
cell with a total of 12 possible vibrational modes. Among these modes only two are Raman active with Eg 
and Ag symmetry. The Ag mode corresponds to vibrations of the Cu-0 bonds along the c-axis while the Eg 
mode represents vibrations in the triangular lattice perpendicular to the c-axis. The atomic displacements 
for these modes are illustrated in Ref . [3S] . In order to determine the symmetry of the modes observed in 
CuFe02 (Fig. [2]) and CuCr02 (Fig. [s]), we performed polarized Raman scattering measurements at room 
temperature using two different laser sources. Here, Raman scattering geometries are identified using the 
Porto notation ki{eies)ks. The labels z' and y' designate directions making an angle 6 relative to the z and 
y axes, where 9 = 50° for CuFe02 while 6 = 15° in the case of CuCr02. 

To our knowledge, no polarized Raman measurements on CuFe02 single crystals have been reported. At 
room temperature, the spectrum taken with the Ar+ laser using unpolarized (u) scattered light, Fig. J2k, 
shows modes at 349 cm~^ and 690 cm~^ in good agreement with results obtained on polycrystals [28l 129) . 
The intensity of the mode at 690 cm^^ disappears with cross (y'x) polarization while the mode at 349 cm~^ 
remains visible in both (y'u) and cross (y'x) polarizations. Measurements with the He-Ne laser show Raman 
modes at 351 cm~ and 692 cm~ and a broad band at 496 cm" -1 (See Fig. Dd). The mode at 692 cm 
has a strong intensity in the parallel polarization (yy) and disappears in the cross polarization (yx). The 
intensity of the mode at 351 cm~^ is very weak which implies that the He-Ne excitation line at 632.8 nm is 
not in resonance with the vibrations associated with this mode as observed in LiNi02 [26] . Despite its weak 
intensity, it is visible in both polarizations. Moreover, this mode was reproducible down to low temperatures 
(see Fig. |4]d). According to the Raman scattering tensors associated with the trigonal point group 3m [57], 
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a cross polarization configuration such as z{yx)z allows only Eg modes, while a parallel polarization config- 
uration like z{xx)z allows the observation of Eg and Ag modes. Therefore, the mode symmetry is assigned 
as WA — 692 cm""'^ and we — 351 cm~^. 

In addition to the vibrational modes observed in CuFe02, a broad band located at 496 cm""'^ is also 
revealed using both laser sources. In the unpolarized {y'u) spectrum obtained with the Ar+ laser, the 
intensity of this feature is within the background noise. In the parallel polarized spectrum obtained with 
the He-Ne laser, the broad peak is clearly observed and detected down to 5 K (Fig. |4|. Polarized spectra 
obtained with both excitation lines show that the mode at 496 cm"^ has an Ag symmetry. The possible 
origin of this band will be discussed later. 

As in the case of CuFe02, CuCr02 should show two Raman modes. However, with unpolarized [u] 
scattered light (not shown) or a parallel polarization [xx) configuration with the Ar"*" laser (Fig. |3^), we 
observe modes at 104 cm^^, 207 cm~^, 382 cm~^ 457 cm~^, 538 cm~^, 557 cm~^, 623 cm~^, 668 cm~^, 
and 709 cm~^. Using the He-Ne laser (Fig. |3]d), a similar spectrum is obtained except for a mode at 382 
cm^^ using the Ar+ laser and a mode at 359 cm~^ in the case of the He-Ne laser. The symmetries of these 
modes can be assigned according to the polarized Raman measurements shown in Fig. [3] Since the modes 
at 104 cm~^, 212 cm~^, and 457 cm~^ are observed in both parallel and cross polarized Raman spectra, 
these modes belong to the Eg irreducible representation (IR) . The other modes are therefore assigned to Ag 
IR since their intensities either are weak or disappear in the cross polarization configuration. So far, there 
have been four publications reporting Raman spectra on CuCr02 powder samples [30], l3T] l32] l33] . Two of 
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Figure 2: Polarized Raman Spectra of CuFe02 at room temperature obtained using the Ar+ and He-Ne 
lasers. Experimental scattering geometries are represented by the Porto notation above each spectrum. 
Polarizations along the y' and x axes arc parallel and perpendicular to the plane of incidence, respectively. 
Strong plasma line at 521 cm~^ in the z'{y'u)z' spectrum was removed for clarity. Raman modes with Kg 
and Eg symmetries are located at = 692 cm~^ and coEg =351 cm~^, respectively. 



5 



these publications show modes at 207 cm~^, 444 cm"^, and 691 cm"^ [301 13I]- In addition, one of these 
works shows additional features with weak intensities at 540 cm"-'^ and 560 cm"-'^ [21]. Other publications 
|32(l33j reveal Raman modes only at 452 cm~^ and 703 cm~^. By comparison, our polarized Raman results 
indicate that the Raman modes in CuCr02 correspond to ujp^ — 709 cm^^ and oje — 457 cm^^. 
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Figure 3: Polarized Raman Spectra of CuCr02 at room temperature obtained using the Ar+ and He-Ne 
lasers. Experimental scattering geometries are designated by the Porton notation above each spectrum. 
Polarizations along the y' and x axes are parallel and perpendicular to the plane of incidence, respectively. 
Green arrows indicate the plasma lines (P). Raman modes have frequencies at we — 457 cm~^ and wa = 
709 cm-i. 

As mentioned earlier, CuFe02 and CuCr02 should only have two Raman modes. However, both com- 
pounds show additional features (Figs. [2] and |3| similar to those observed in other delafossite compounds 
such as CuA102 ^34j and CuGa02 |25j . In agreement with ab initio calculations, these additional modes in 
CuA102 are attributed to non-zero wavevector phonons which are normally forbidden by Raman selection 
rules [31]. As suggested, the selection rules are possibly relaxed by defects such as Cu vacancies, interstitial 
oxygens or tetrahedrally coordinated Cr+^ or Fe"'"'^ on the Cu site [Mj. Thus, the additional features ob- 
served in CuFe02 and CuCr02 could have an origin similar to that observed in CuA102 [Ml and CuGa02 
[25j . They could also be related to crystal field excitations which were revealed in Raman spectra of other 
geometrically frustrated compounds |35j . 

3.2 Temperature dependent measurements 

Unpolarized Raman spectra of CuFe02 obtained between 290 K and 5 K are presented in Fig. |4^. Over this 
temperature range, no considerable change is observed. In particular, no splitting of the Eg mode below T^i 
is noticeable despite the R3m C2/m structural transition at T/vi [I1[71|T7]. We attribute this discrepancy 
to weak resonance with the He-Ne excitation line, which results in the weak intensity of the mode and 
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makes it difficult to resolve any possible splitting. Another possibility is that the temperature of the sample 
remains above T^i even with a beam power of 4000 W/cm^. 




Figure 4: Raman spectra of (a) CuCr02 and (b) CuFe02 between 5 K and 290 K. While no considerable 
change is observed in the CuFe02 spectra in this temperature range, an additional mode in CuCr02 appears 
at 467 cm~^ below 200 K and slightly increases in frequency at lower temperatures. In addition, Raman 
spectrum of CuCr02 at 8 K shows a broad band centered at 550 cm"-'^. This mode could be due to magnon 
modes reflecting the proper screw spin structure below T/v2 = 23.6 K in CuCr02. 



In the case of CuCr02, unpolarized Raman spectra shown in Fig. display noticeable differences as 
the temperature is decreased from room temperature down to 8 K. With a close look at the Eg mode at 
458 cm^^, one can observe that its tale becomes broader on the right hand side starting at 200 K. With 
further cooling, an additional mode is easily distinguished and its frequency increases to 470 cm^^ at 8 K. 
This mode is also observed with parallel and cross polarization configurations (not shown). It should be 
noted that neutron diffraction measurements [3 117j , magnetostriction measurements [23 and sound velocity 
measurements (Fig. [T]) do not show any anomaly that could be associated with a structural deformation in 
the temperature range around 200 K. This mode could have an origin similar to that of the additional modes 
observed between room temperature and 8 K (Fig.[4]D). Moreover, the spectrum of CuCr02 at 8 K (Fig. |4]d) 
deserves some attention. Unlike the spectra at other temperatures, it develops a broad background feature 
centered at ^550 cm"-'^. This broad band, which can also be observed using parallel and cross polarizations 
(not shown), might be associated with magnon modes owed to a proper screw ordering observed below Tn2 
[36II33. Finally, although there is some evidence for a structural deformation at TVi in CuCr02 [23], no 
additional Raman modes are observed below this temperature. Although local heating due to incident beam 
power (4000 W/cm^) is possible, the broad band observed in the spectrum at 8 K (Fig. j^ja)) clearly shows 
that the temperature is below T/vi. Linewidths of Raman modes normally narrow down with decreasing 
temperature. For example, in BiFeOa, which undergoes a structural transition at the Curie temperature Tc 
= 1100 K, Fukumura et. al. |38j observed only 7 Raman modes at room temperature due to broadening of 
the modes. All 13 Raman modes were observed only at 4 K much below [35] • For CuCr02, even lower 
temperatures and a lower incident beam power might be required for the observation of additional modes. 
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Temperature variations of the frequencies of the Raman modes in CuFe02 and CuCr02 are presented 
in Fig. [5] As shown in Fig. [5^, the frequencies of both modes in CuFe02 increase ahiiost hnearly down to 
50 K with no significant variation below T/vi — 14 K. Similarly, in the case of CuCr02 (Fig. Isb), the mode 



frequencies increase between 290 K and 80 K. While the Ag mode frequency remains constant between 80 
K and 8 K, the frequency of the Eg mode seems to decrease slightly in this range. The additional modes 
(Fig. |4| observed in both compounds behave similarly with temperature (not shown). 
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Figure 5: Temperature dependencies of the frequencies of the Raman modes in a) CuFe02 and b) CuCr02 
between 295 K and 5 K could be associated with lattice contraction and anharmonic phonon-phonon inter- 
actions. Therefore, the Eg modes in these compounds are not associated with the order parameters of the 
pseudoproper ferroelastic transitions at T/vi in CuFe02 and CuCr02. 



As discussed earlier, neutron [T7] and x-ray [7] diffraction measurements on CuFe02 reveal an R3m 
C2/m structural transition while magnetostriction measurements on CuCr02 show evidence for crystal 
symmetry lowering. In accordance with these results, sound velocity measurements on CuFe02 [I] and 
CuCr02 (see Fig. [l]) indicate that both compounds undergo an i?3m C2/m pseudoproper ferroelastic 
transition at Tni. According to the group theory [T5], one possible scenario is that an Eg-symmetric optic 
mode is associated with the order parameter [T] . In this case, the excess Gibbs free energy Ge can be written 

as 111 

Ge = -mwo'^u^ + -bu^ + -Gel + jBsU, (3) 

where 

mwo^ = a{T - To) A. (4) 

In the above equations, m is the reduced mass, Wo is the uncoupled frequency of the soft Eg mode, while a 
and b are temperature independent constants. The first two terms in Eq. [3]are due to the Landau expansion 
of the order parameter w, which corresponds to normal coordinate vibrations associated with the soft Eg 
optic mode. Thus, the first term ^mwo^u^ corresponds to the harmonic oscillator energy and is the only 
temperature dependent term. The term ^Ce^ is the elastic energy of the soft acoustic mode associated 
with the strain component eg. For simplification, only one elastic constant G is considered (See Ref. pi for 
complete elastic energy). Finally, the bilinear coupling term jegU in Eq. ^ with 7 representing the coupling 
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coefficient, is necessary in order to account for the softening of the acoustic modes observed in CuFe02 [I] 
and CuCr02 (see Fig. [T]). Minimizing Gg with respect to u and Cs, one obtains 

e« ~ 

C 

and 



(5) 



— , (6) 



showing that the biUnear couphng term {"jegu) renormahzes the uncoupled transition temperature to 

Tni=To + ^. (7) 
Finahy, the frequency uj of the soft optical mode can be obtained using [39] 

which yields 

2 

uj^ = -iT-TNi) + ^ (T > Te) (9) 
m Cm 

and 

io' = —{T-T^,) + ^ (T < T,). (10) 
m Cm 

According to Eqs. [9] and [Toj the frequency square of the soft optical mode should vary linearly with 
temperature with a slope change at T^vi as observed in some pseudoproper ferroelastic compounds |19[ 
[201 [3TJ [22]. According to our Raman measurements, the temperature dependence of the Eg symmetry 
modes cannot be associated with that of a soft optic mode. Thus, the temperature behavior of all modes is 
rather attributed to thermal contraction and anharmonic phonon-phonon interactions, in agreement with the 
analyses of Pavunny et. al. [29] for their Raman study of CuFe02 between 400 K and 80 K. Our conclusion 
is that none of the Raman modes in CuFe02 and CuCr02 can account for the pseudoproper ferroelastic 
transitions observed at T/vi = 14 K in CuFe02 and at T/vi = 24.3 K in CuCr02. While these results do 
not refute the pseudoproper ferroelastic transitions in CuFe02 [Ij and CuCr02 (see Fig. [T]), they leave the 
driving mechanisms unresolved. 



4 Conclusions 

Polarized Raman scattering measurements were performed on delafossite magnetoelectric CuFe02 in order 
to determine the true nature of the order parameter associated with pseudoproper ferroelastic transition 
observed in CuFe02 by means of sound velocity measurements [1]. As preliminary sound velocity mea- 
surements on the isostructural compound CuCr02 show similar elastic softening at the antiferromagnetic 
transition (see Fig.[T]), the Raman measurements were also performed on CuCr02 for comparison. 

Apart from the vibrational modes in CuFe02 and CuCr02 with A^ and Eg symmetries, one additional 
mode in CuFe02 and seven additional modes in CuCr02 are observed at room temperature. Below 200 K, 
another mode in CuCr02 with a frequency close to that of the Eg mode appears and persists to lower tem- 
peratures. The additional modes observed in both compounds are possibly associated with either relaxation 
of Raman selection rules or crystal field excitations. More interestingly, the spectrum of CuCr02 at 8 K 
shows a broad band centered at 550 cm"-'^ attributed to the proper screw ordering below Tpf2 in CuCr02. 
Furthermore, for both compounds the frequency of all modes increase with decreasing temperature, with 
no significant variation at the phase transitions. The observed temperature behavior is thus attributed to 
thermal contraction and anharmonic phonon-phonon interactions. Therefore, these results show that the Eg 
symmetry Raman modes in CuFe02 and CuCr02 do not induce the transitions observed at Tjvi, leading to 
the necessity of further search for the true origin of the order parameters associated with the pseudoproper 
ferroelastic transitions observed at Tjvi in both compounds. 
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